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Abstract. The operation of binary intermolecular recombination, originating 
in the theory of DNA computing, permits a natural generalization to ra-ary 
operations which perform simultaneous recombination of n molecules. In the 
case n = 3, we use computer algebra to determine the polynomial identities of 
degree < 9 satisfied by this trilinear nonassociative operation. Our approach 
requires computing a basis for the nullspace of a large integer matrix, and for 
this we compare two methods: (i) the row canonical form, and (it) the Hermite 
normal form with lattice basis reduction. In the conclusion, we formulate some 
conjectures for the general case of n-ary intermolecular recombination. 



1. Introduction 

Let X be a finite set, and let X* be the free semigroup on X; thus X* consists 
of all finite nonempty strings of elements of X, 

X* = { XiX 2 •■•Xk\xi,X2,...,Xk <E X (k > 1) }, 

together with the binary operation of concatenation, 

(Xl#2 • ' • Xk){Xk+lX k +2 ■ ■ ■ Xk+e) 1 — > xix 2 ■ ■ ■ Xk+e- 

Let n be a positive integer, and let (X*) n be the set of all n-tuples of elements of 
X*. Let R(X, n) be the vector space with basis (X*) n over the field Q of rational 
numbers. 

Definition 1. On R(X,n) we introduce a multilinear nonassociative operation 
{a\, 0,2, • - • ) a n }> called n-ary intermolecular recombination, and defined on 
basis elements as follows: For = (an, a 2 2, . . . , a,i n ) G (X*) n (i = 1, . . . , n) we set 

(1) {ai,a 2 , . . . ,a„} = ^ (<V(1),1, 0<r(2),2) • • • ? «a-(n),n)) 

<res„ 

where S n is the symmetric group on {1, 2, . . . , n}. 

In Definition [1] each at can be regarded as the abstract representation of a 
molecule divided into n consecutive submolecules an, a i2 , . . ., ai n . For example, if 
X = {A, C, G, T} represents the four bases found in DNA, then each a^ represents a 
DNA sequence partitioned into n subsequences. With this interpretation, the sum 
defining {ai, a 2 , . . . ,a n } expresses the results of recombining the n molecules in all 
possible ways which preserve the positions of the n submolecules in each molecule. 
Formula ([1} is very similar to the definition of the permanent of an n x n matrix; 
see Wanless [?]■ 

The next result is clear. 
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Lemma 2. The operation of n-ary intermolecular recombination in Definition^ 
satisfies the complete symmetry identity in degree n: 

(2) {<V(1)> a 7r(2)j • • ■ j = W) 0,2i ■ • ■ 7 0>n}, 

for every tt G S^; i/ins {ai,ci2, . . . ,a n } depends only on the multiset of arguments, 
not the order of the arguments. 

In the case n = 2, the operation of binary intermolecular recombination was 
introduced by Landweber and Kari [J. 

Bremner [I] used computer algebra to establish the following proposition. 

Proposition 3. Every polynomial identity of degree < 5 satisfied by binary inter- 
molecular recombination is a consequence of commutativity (binary complete sym- 
metry) and the binary recombination identity 

2{{{a, b}, c}, d} - {{{a, b}, d}, c] - {{{a, c}, b}, d} - {{{6, c}, a}, d} 

+ {{a,b},{c,d}} = 0. 

Sverchkov [5] has recently proved the following theorem. 

Theorem 4. Every identity, with no restriction on the degree, satisfied by binary 
intermolecular recombination is a consequence of commutativity and the binary re- 
combination identity. 

In related work, Bremner, Piao and Richards 3 have shown that every polyno- 
mial identity satisfied by the zygotic algebras of simple Mcndelian inheritance is a 
consequence of commutativity and the binary recombination identity. 

In the present paper we consider the case n = 3; we use computer algebra to de- 
termine the polynomial identities of degree < 9 satisfied by ternary intermolecular 
recombination. A monomial which involves d applications of an n-ary operation 
has degree d(n— 1) + 1; thus polynomial identities for ternary intermolecular recom- 
bination exist only in odd degrees. 

2. Binary intermolecular recombination 

In this section we recall the results of Bremner [1] on polynomial identities for 
binary intermolecular recombination, and show how these results may also be ob- 
tained (and slightly improved) using the Hermite normal form with lattice basis 
reduction. 

In degree 4, there are 15 commutative nonassociative multilinear monomials, 12 
with association type {{{a, &}, c}, d} and 3 with association type {{a, &}, {c, d}}. 
We order these monomials first by association type, and then lexicographically by 
the permutation of the variables: 

{{{a,b},c},d}, {{{a,b},d},c}, {{{a,c},b},d}, {{{a, c}, d}, b}, 

{{{a,d},b},c}, {{{a,d},c},b}, {{{b,c},a},d}, {{{6, c}, d}, a}, 

{{{b,d},a},c}, {{{b, d},c}, a}, {{{c,d},a},b}, {{{c, d}, b}, a}, 

{{a,b},{c,d}}, {{a,c},{b,d}}, {{a,d},{b,c}}. 

We expand each monomial by three applications of binary intermolecular recombi- 
nation. For example, writing 

a = (ax,a 2 ), b=(b 1 ,b 2 ), c=(ci,c 2 ), d=(di,d 2 ), 
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Table 1 . The expansion matrix E in degree 4 for n = 2 



we obtain 

{{{a,b},c},d} = (ai,d 2 ) + (&i,d 2 )+2(ci,d 2 ) + (di,a 2 ) + (d x , b 2 ) + 2(di, c 2 ), 
{{a, b}, {c, d}} = (ai,c 2 ) + (ai,d 2 ) + (61,02) + (&i,d 2 ) + (ci,a 2 ) + (ci,6 2 ) 
+ (di,a 2 ) + (di,b 2 ). 

Each expansion is a linear combination (allowing zero coefficients) of 12 pairs, which 
we order lexicographically: 

(ai,6 2 ), (ai,c 2 ), (ai,d 2 ), (61,02), (61, c 2 ), (61,^2), 

(ci,a 2 ), (ci,6 2 ), (ci,d 2 ), (di,a 2 ), (^1,62), (di,c 2 ). 

We store the expansions in the 12 x 15 matrix E in which entry contains the 
coefficient of pair i in the expansion of monomial j; see Table [TJ The coefficient 
vectors of the polynomial identities satisfied by binary intermolecular recombination 
in degree 4 are the nontrivial linear combinations of the columns of E that give the 
zero vector; that is, the nonzero vectors in the nullspace of E. 

There are two ways to compute a basis for the nullspace of a matrix E with 
entries in the ring Z of integers: 

(a) Compute the row canonical form (RCF) of E over Q, set the free variables 
equal to the standard basis vectors in Q d where d is the dimension of the 
nullspace, and solve for the leading variables; then (if necessary) multiply 
each vector by the LCM of the denominators of its components to obtain 
integer components and divide each new vector by the GCD of its (now 
integral) components. 

(b) Compute the Hermite normal form H of the transpose E l together with a 
transform matrix U for which UE 1 = H. The bottom d rows of U form a 
basis for the integer nullspace lattice of E. Then use the LLL algorithm for 
lattice basis reduction on the bottom d rows of U. (This method was first 
applied to polynomial identities by Bremner and Peresi [2].) 

For the matrix E of Table [TJ method (a) gives the following results. The RCF 
over Q is given in Table [5J and the corresponding nullspace basis vectors (sorted 
by increasing length) are given in Table [3] The squared Euclidean norms of these 
basis vectors are 6, 8, 8, 12, 20, 22, 34, 52, 70. 
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Table 2. The row canonical form of the matrix E for n = 2 
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Table 3. The canonical basis of the nullspace of E for n = 2 



For each basis vector, we apply all permutations of a, b, c, d to the corresponding 
identity / and store the results in a 24 x 15 matrix M(J); the rank of M(I) is the 
dimension of the ^-module of identities which are consequences of I. Rows 2 and 
3 are the shortest vectors (squared norm 8) for which the corresponding identities 
generate the entire nullspace, and row 3 corresponds to the binary recombination 
identity. 

Method (a) provides a basis of integral vectors for the rational nullspace of the 
matrix E. The disadvantage of this method is that this basis of integral vectors 
may not be a integer basis for the nullspace lattice of E] see Example 17 in Bremner 
and Peresi [2J. 

Definition 5. The nullspace lattice of the m x n integer matrix E is 

L(E) = {X e Z" I EX = O}. 

We want to find a lattice basis of L(E); that is, a set of vectors in L(E) which 
are linearly independent over Q and such that every vector in L(E) is an integer 
linear combination of these basis vectors. To find a lattice basis, we need to recall 
the definition of the Hermite normal form (HNF) of an integer matrix. 

Definition 6. The mxn matrix H over Z is in Hermite normal form (HNF) if 
there exists an integer r (the rank of H) with < r < m and a sequence of integers 
1 < ji < j% < • • • < ir < n such that 

(1) Hij = for 1 < i < r and 1 < j < j it 

(2) H iJz > 1 for 1 < i < r, 

(3) < /// . , < H idi for 1 < i < r and 1 < k < i, 

(4) : = for r + 1 < i < m and 1 < j < n. 

The following two results are Lemmas 19 and 20 in Bremner and Peresi [2J. 
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Table 4. Hermite normal form H of E l for n = 2 
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Table 5. Transform matrix U with UE t — H for n = 2 



Lemma 7. If E is an m x n matrix over Z, then there is a unique m x n matrix 
H over Z in HNF such that UE = H for some m x m matrix U over Z with 
det(U) — ±1. (The matrix U is not unique.) 

Lemma 8. Let E be an m x n matrix over 1, let H be the HNF of E l , and let U 
be an n x n matrix over Z with det(C/) = ±1 and UE 1 — H. If r is the rank of H , 
then the last n — r rows of U form a lattice basis for L(E). 

For the matrix E of Table [1] method (b) implemented with the Maple command 

HermiteForm( Transpose (E) , output=['H' , 'U'] ): 

gives the Hermite normal form H in Table 0] and the transform matrix U in Table 
O Since the rank of E is 6, the last 9 rows of U form a lattice basis of the integer 
nullspace of E l . The squared norms of these 9 row vectors are 11, 5, 5, 31, 15, 25, 
6, 124, 44, 196, 14, 142, 22, 6, 38. At this point the lengths of the basis vectors 
are greater than those obtained with the RCF; however, we know that we have an 
integer basis of the nullspace lattice. 
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Table 6. Reduced transform matrix U with UE l = H for n = 2 



To improve these results we use the LLL algorithm for lattice basis reduction, 
following Brcmncr and Peresi [2]. For the matrix E of Table [1] the Maple command 

HermiteForm( Transpose (E) , output='U' , method= ' integer [reduced] ' ): 

produces the transform matrix U of Tabled The bottom 9 rows of U are a reduced 
basis for the integral nullspace lattice of E; every one of these vectors has squared 
norm 6, the same as the shortest vector from Table [3] Even for this small matrix, 
the LLL algorithm has produced a basis of the nullspace with significantly shorter 
vectors. The identity corresponding to row 9 of U generates the entire nullspace: 

{{{a, b}, c}, d} - {{{a, c}, b}, d} - {{{6, c}, d}, a} - {{{6, d}, a}, c} 
+ {{{b,d},c},a} + {{a,b},{c,d}} = 0. 

Identity (j4]) has one more term than the binary recombination identity |3j), but its 
coefficient vector is shorter. 

Identities and ^ both have only one term in the second association type, 
so both imply that every monomial in the second type can be expressed as a linear 
combination of monomials in the first type. For example, identity (j4]) can be written 
in the form 

{{a, b}, {c, d}} = - {{{a, b}, c}, d} + {{{a, c}, 6}, d} + {{{6, c}, d}, a} 
+ {{{b,d} 1 a},c}-{{{b,d} 1 c},a}. 

This identity can be used to reduce the number of association types that we need 
to consider when studying identities of higher degree. 

3. Ternary intermolecular recombination 

We now present the main results of this paper: a complete and minimal set of 
polynomial identities of degree < 9 for ternary intermolecular recombination. 
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3.1. Degree 3. Every polynomial identity in degree 3 satisfied by ternary inter- 
molecular recombination is a consequence of the complete symmetry identity. In- 
deed, equation ^ implies that there is only one monomial in degree 3, namely 
{a, b, c}, and hence the only possibly polynomial identity in degree 3 is {a, 6, c} = 
for all a, b, c; this clearly does not hold. 

3.2. Degree 5. In degree 5 there are three distinct association types for a ternary 
operation: 

{{a,b,c},d,e}, {a, {b, c, d}, e}, {a, b, {c, d, e}}. 
Equation ^ implies that we need to consider only the first type, and that we 
need to consider only 10 distinct monomials in this type; we order these monomials 
lexicographically : 

{{a, b, c},d, e}, {{a, b, d}, c, e}, {{a, 6, e}, c, d}, {{a, c, d}, b, e}, 

{{a,c,e},6,e}, {{a, d, e}, b, c}, {{b, c, d}, a, e}, {{6, c, e}, a, d}, 

{{&, d, e},a, c}, {{c,d,e},a, 6}. 

Any polynomial identity in degree 5 satisfied by ternary intermolecular recombina- 
tion is a linear combination of these 10 monomials. 

Lemma 9. Every polynomial identity of degree 5 satisfied by ternary intermolecular 
recombination is a consequence of the complete symmetry identity in degree 3. 

Proof. Consider the expansions of the 10 monomials in degree 5 using two applica- 
tions of ternary intermolecular recombination. For example, since 

{a, 6, c} = (ai,6 2 ,c 3 ) + (ai,c 2 ,6 3 ) + (&i,a 2 ,c 3 ) 

+ (6i,C2,o 3 ) + (01,02,63) + (ci,6 2 ,a 3 ), 

we obtain 

{{a, 6, c}, d, e} = 2(ai, d 2 , e 3 ) + 2(ai, e 2 , 4) + 2(di, a 2 , e 3 ) + 2(d x , e 2 , 03) 

+ 2(ei,a 2 ,d 3 ) + 2(ei,d 2 ,a 3 ) + 2(6i,d 2 ,e 3 ) + 2(6i,e 2 ,d 3 ) 

+ 2(d 1 ,b 2) e 3 ) + 2(d u e 2 , 63) + 2(e x , 6 2 , d 3 ) + 2( ei ,d 2 ,b 3 ) 

+ 2(ci,d 2 , e 3 ) + 2(ci,e 2 ,d 3 ) + 2(di,c 2 ,e 3 ) + 2(di,e 2 ,c 3 ) 

+ 2(ei,c 2 ,d 3 ) + 2(ei,d 2 ,c 3 ). 

Each term in these 10 expansions has the form (%i,y 2 , z 3 ) where (x,y, z) is a per- 
mutation of a 3-element subset of {a, 6, c, d, e}. There are 60 such triples, which we 
order lexicographically: 

(ai,6 2 ,c 3 ), (ai,6 2 ,d 3 ), (ai,6 2 ,e 3 ), (a x ,c 2 , 63), (ai,c 2 ,d 3 ), (ai,c 2 ,e 3 ), 

(ei,c 2 ,a 3 ), (ei,c 2 ,6 3 ), (ei,c 2 ,d 3 ), (ei,d 2 ,a 3 ), (ei,d 2 ,6 3 ), (ei,d 2 ,c 3 ). 

We construct the 60 x 10 matrix E in which entry contains the coefficient 

of triple i in the expansion of monomial j. The polynomial identities in degree 5 
for ternary intermolecular recombination correspond to the nonzero vectors in the 
nullspace of E. It suffices to show that this matrix has full rank, and for this it 
suffices to find 10 rows for which the corresponding 10 x 10 submatrix has full rank. 
Rows 1, 2, 3, 5, 6, 9, 17, 18, 21 and 33 produce the submatrix displayed in Table 
the row canonical form of this submatrix is the identity matrix. □ 
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Table 7. Submatrix for the proof of Lemma |9] 



3.3. Degree 7. We now consider 7 distinct molecules, each divided into 3 sub- 
molecules: 

a = (at, 0,2,0,3), b = (61,62,63), c = (ci,c 2 ,c 3 ), d = (d 1} d 2 , d 3 ), 

e = (ei,e 2 ,e 3 ), / = (/ x , / 2 , / 3 ), 9^(91,92,93)- 

After performing three applications of ternary intermolecular recombination on 
these molecules with some order of operations and some permutation of molecules, 
we obtain a linear combination of triples of the form (x±, y 2 , z 3 ) where (x, y, z) is a 
permutation of a 3-element subset of {a, b, c, d, e, /, 5}. The number of triples which 
can occur as terms in these linear combinations is therefore 

81 Q =210. 

We order these triples lexicographically by the permutation (x, y, z). 

Equation ([2]) implies that we need only two association types in degree 7, 

{{{a,b,c},d,e}J,g} and {{a, b, c}, {d, e, /}, g}, 

which contain respectively the following numbers of distinct monomials, 

7! 1 7! 

= 210 and — • — — = 70. 



3!2!2! 2! 3!3!1! 

Within each association type, we order the monomials lexicographically by the per- 
mutation of the variables. Any polynomial identity in degree 7 satisfied by ternary 
intermolecular recombination is a linear combination of these 280 nonassociativc 
monomials. 

The expansion matrix E in degree 7 has 210 rows and 280 columns; entry 
is the coefficient of triple i in the expansion of nonassociative monomial j. The 
expansions of the two association types are displayed in Tables [8] and [9] Each 
expansion has 6 3 = 216 terms; after collecting terms, the first type produces a 
linear combination of 30 triples with multiplicities 4 and 12, and the second type 
produces a linear combination of 54 triples with multiplicity 4. We use Maple to 
find that the matrix E has rank 35, and so its nullspace has dimension 245. We 
compute the row canonical form of E, and find that every entry of the RCF is 
an integer. We obtain the canonical basis of the nullspace of E; this is a list of 
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4[ai,/ 2 ,53]+ 4[oi,sa J /s]+ 4[6 X , / 2 , 53] + 4[&i,ft,,/ 3 ] + 4[ Cl ,/ 2 , 33 ] 
+ 4[ci,sa, / 3 ] + 12 [di, / 2 , 9s] + 12[di,fla, / 3 ] + 12[ei, /2,9a] + 12 [ei, 92, /s] 
+ 4[/i, a 2 ,9 3 ]+ 4[/i,& 2 ,93] + 4[/i,c 2 ,9 3 ] +12[/i,d 2 ,93] + 12[/i, e 2 , 93] 
+ 4[/i, 92, 03] + 4[/i, ff2 , 63] + 4[/i, g 2 , cs] + 12[/i, ff2 , <fe] + 12[/i, 92, e 3 ] 
+ 4[9i, a 2 , f 3 ] + 4[9i, 62, /s] + 4[ 3l , c 2 , / 3 ] + 12[.gi, d 2 , / 3 ] + 12[ 5 i, e 2 , / 3 ] 
+ 4[3i,/ 2 ,a 3 ]+ 4[9i,/ 2 ,6 3 ]+ 4[pi, / 2 , c 3 ] + 12[.gi, f 2 , d 3 ] + 12 [91, / 2 , e 3 ] 
Table 8. Expansion of monomial {{{a, 6, c}, d, e}, /, 3} 

4[ai,d 2 ,3 3 ] + 4[ai,e 2 ,3 3 ] +4[ai,/ 2 ,3 3 ] +4[ai,3 2 ,d 3 ] +4[ai,3 2 ,e 3 ] 
+ 4[ai,3 2 ,/ 3 ] + 4[& 1 ,d 2 ,9 3 ] + 4[&i,e 2 ,9 3 ] + 4[& 1 ,/ 2 ,9 3 ] + 4[6i,9 2 ,d 3 ] 
+ 4[6i,92,e 3 ] + 4[6i,92,/ 3 ] + 4[ci,d2,9 3 ] + 4[ci, e 2 , 33] + 4[ci, / 2 , 33] 
+ 4[ci,92,d 3 ] + 4[ci,3 2 ,e 3 ] + 4[c x , 92, / 3 ] +4[di, a 2 , 93] + 4[di, 6 2 , 93] 
+ 4[di,c 2 ,3 3 ] +4[di,92,a 3 ] +4[di,9 2 ,& 3 ] + 4[di, 92, c 3 ] +4[e 1 ,a 2 ,9 3 ] 
+ 4[ei,6 2 ,3 3 ] + 4[ei,c 2 ,3 3 ] +4[ei,3 2 ,a 3 ] + 4[ei,9 2 ,6 3 ] + 4[ei,9 2 ,c 3 ] 
+ 4[/i,a 2 ,3 3 ] + 4[/i,6 2 ,3 3 ] + 4[/i, c 2 , 33] + 4[/i, 32, a 3 ] + 4[/i,3 2 ,6 3 ] 
+ 4[/i, 92, c 3 ]+4[3x,a2,d 3 ] + 4[9i,a 2 ,e 3 ]+4[3x,a2,/3] + 4(91,62, 4] 
+ 4[3i,6 2 ,e 3 ] + 4 [31,62,73] +4[9i,c 2 ,d 3 ] + 4 [31, c 2 , e 3 ] + 4 [91, c 2 , / 3 ] 
+ 4[3i,d 2 ,a 3 ] + 4[3i,d 2 ,6 3 ] +4[3i,d 2 ,c 3 ] +4[9i,e 2 ,a 3 ] + 4[9i,e 2 ,6 3 ] 
+ 4[9i,e 2 ,c 3 ] +4[9i,/ 2 ,a 3 ] + 4[3i,/ 2 ,6 3 ] + 4[3i,/ 2 ,c 3 ] 

Table 9. Expansion of monomial {{a, 6, c}, {d, e, /}, 3} 



245 vectors of dimension 280 with integer components. We sort these vectors by 
increasing Euclidean norm; the list of squared norms is displayed in Table 1101 

We now perform further computations, using modular arithmetic to save mem- 
ory, to determine a set of generators for the nullspace as a module over the sym- 
metric group 6*7. (We use p = 101, the smallest prime greater than 100; any prime 
greater than the degree of the identities would produce the same dimensions.) For 
each basis vector in the nullspace, we apply all 5040 permutations of a, 6, c, d, e, /, 
9 to the corresponding identity, and store the results in a matrix of size 5040 x 280. 
We compute the row canonical form of this matrix; its nonzero rows form a basis 
of the SV-module generated by the identity, and its rank is the dimension of this 
module. We process the nullspace basis vectors in order, saving previous results 
(the nonzero rows which form a basis of the SV-module generated by the previous 
basis vectors). At each stage, an increase in the rank implies that the current iden- 
tity is a new generator of the nullspace as a module over S7. The results of these 
computations show that the canonical basis vectors with positions 1, 10 and 60 
(after the basis vectors have been sorted by increasing Euclidean norm) represent 
polynomial identities which generate the entire nullspace as an SV-module: 



1 = {{{a,6,c},d,e},/,3} - {{{a, 6, e}, c, d}, /, 9} - {{{6, c,d} 7 a,e}, f, g} 
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1386 1386 1402 1402 1408 1412 1414 1414 1440 1470 

1482 1484 1494 1516 1568 1592 1678 1758 1818 1824 

1852 1894 1908 1910 1918 1946 1946 1982 2002 2026 

2064 2174 2196 2226 2246 2272 2312 2324 2326 2340 

2346 2366 2390 2396 2402 2418 2430 2472 2490 2562 

2634 2640 2654 2668 2844 

Table 10. Squared norms of canonical basis vectors in degree 7 
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Table 1 1 . Squared norms of reduced basis vectors in degree 7 
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+ {{{b,d,e},a 7 c},f,g}, 
J = {{{a,b,e},c,g},d, /} - {{{a, b, e}, /, g}, c, d}- {{{a, c, e}, b, g}, d, /} 

+ {{{a, c, e}, /, 5}, 6, d} + {{{a, e, /}, 6, 5}, c, d} - {{{a, e, /}, c, 5}, 6, d}, 
K = 2{{{a, b, c}, d, e}, f, g} - {{{a, b, c}, d, g}, e, /} - 2{{{a, 6, c}, e, 5}, d, /} 

+ 2{{{a, 6, d}, c, e}, /, g} - 2{{{a, b, d}, c, g}, e, /} - {{{a, 6, d}, e, <?}, c, /} 

- {{{a, b, e}, c, d}, f, g} + {{{a, b, e}, c, g}, d, /} + 3{{{a, 6, 5}, c, d}, e, /} 

- {{{«, c, d}, 6, e}, /, 3} - {{{6, c, d}, a, e}, /, 5} + {{a, fc, e}, {c, d, g}, /}. 

Every identity in the nullspace is a linear combination of permutations of these 
three identities, which have squared norms 4, 6 and 32; the second identity has 
been normalized so that it begins with a positive coefficient. 

We can get significantly better results using the Hermite normal form and lattice 
basis reduction. 

Theorem 10. Every polynomial identity of degree 7 satisfied by ternary inter- 
molecular recombination is a consequence of complete symmetry in degree 3 and 
these three identities in degree 7: 

P = {{{a,c,f},e,g},b,d}~ {{{a, c, g}, e, /}, b, d} - {{{c, e, /}, a, g}, b, d} 

+ {{{ c ,e,g},a, f},b, d}, 
Q = {{{a,d,e},b,g},c, /} - {{{a, d, e}, /, g}, b, c} - {{{b, d, e}, a, g}, c, /} 

+ {{{b, d, e}, /, 5}, a, c} + {{{d, e, /}, a, g}, b, c} - {{{d, e, /}, b, g}, a, c}, 
R = {{{a, b, c}, f, g}, d, e} + {{{a, b, g}, d, /}, c, e} + {{{a, c, d}, b, /}, e, 3} 

- {{{a, c, d}, b, g}, e, /} + {{{a, c, d}, /, g}, b, e} - {{{a, c, /}, b, d}, e, 3} 
+ {{{«, c, g}, d, /}, 6, e} - {{{a, d, /}, c, g}, b, e} - {{{a, d, g}, b, c}, e, /} 

- {{{aJ,g},b, c},d,e} + {{{b,c,d},a,g},e,f} - {{a,d,g},{b,c, f},e}. 

Identities P and Q are independent: neither implies the other. Identity R implies 
both P and Q, but identities P and Q together do not imply R. 

Proof. We apply the Maple command 

HermiteForm( Transpose (E) , output='U' , method= ' integer [reduced] ' ): 

to the transpose (size 280 x 210) of the expansion matrix E, and obtain a transform 
matrix U (size 280 x 280). The bottom 245 rows of U are a basis for the integral 
nullspace lattice of E. We sort these vectors by increasing Euclidean norm; the 
squared norms are displayed in Table 111! It is clear by comparing Tables [10] and 
ITTI that the basis vectors obtained in this way are generally much shorter than 
those obtained from the RCF. Identities P, Q and R have coefficients ±1 and the 
squared norms of the corresponding coefficient vectors are 4, 6 and 12. (These 
identities correspond to the reduced basis vectors in positions 1, 40 and 129 after 
the reduced basis vectors have been sorted by increasing Euclidean norm.) These 
three identities can be verified by expanding each term using three applications of 
ternary intermolecular recombination, and then checking that the results collapse 
to zero. 

To prove the stated dependence and independence relations we proceed as fol- 
lows. Given an identity I (a, b, c, d, e, /, g) we apply all 5040 permutations of a, b, 
c, d, e, /, g and store the coefficients of the resulting identities in a 5040 x 280 
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matrix. The rank of this matrix is the dimension of the SV-submodulc of identities 
generated by /. To check whether another identity J(a, 6, c, d, e, /, g) is implied 
by / we simply determine whether J is in the row space of this matrix. We find 
that identity P produces dimension 105, and identity Q produces dimension 127; 
stacking the two corresponding matrices together, we see that these two identities 
together produce dimension 155. However, identity R produces dimension 245, and 
this is also the rank produced by the three identities together. Since 245 is also the 
dimension of the nullspace of the expansion matrix, it follows that every identity 
in degree 7 is a consequence of the identity R, in the sense that every identity in 
degree 7 is a linear combination of permutations of i?(a, b, c, d, e, /, g). □ 

The identity R of Theorem ITUl will be called the ternary recombination iden- 
tity. The last term in R is the only term in identities P, Q and R which has the 
second association type, so we have the following result. 

Corollary 11. Every polynomial identity of degree < 7 satisfied by ternary inter- 
molecular recombination is a consequence of complete symmetry and the following 
identity, which expresses any monomial in the second association type as a linear 
combination of monomials in the first association type: 

{{a,b, c}, {d,e,f},g} = 

{{{a, d, e}, f, c}, b, g} + {{{a, d, c}, b, /}, e, g} + {{{a, e, b}, d, /}, g, c} 

- {{{a, e, b}, d, c},g, /} + {{{a, e, &}, /, c}, d, g} - {{{a, e, /}, d, b}, g, c} 
+ {{{a, e, c}, b, /}, d, g} - {{{a, b, /}, e, c}, d, g} - {{{a, b, c}, d, e}, g, /} 

- {{{a,f,c},d,e},b,g} + {{{d,e,b},a,c},g, f}. 

(We have permuted the variables so that the term in the second association type has 
the identity permutation.) 

3.4. Degree 9. The number of permutations of 3-element subsets of a set of 9 
variables is clearly 




= 504. 



This is the number of rows in the expansion matrix in degree 9. 

Equation (|2|) implies that we need only four association types in degree 9: 

9' 

{{{{a, b, c}, d, e}, /, g}, h, i} with T^j^y = 7560 monomials, 

1 9! 

{{{a, 6, c}, {d, e, /}, g}, h, i} with — • = 2520 monomials, 

9! 

{{{a, 6, c}, d, e}, {/, g, h}, i} with ^7^7 = 5040 monomials, and 

1 9! 

{{a, b, c}, {d, e, /}, {g, h, i}} with — • = 280 monomials. 

The total is 15400; this is the number of columns in the expansion matrix. 

We initialize the expansion matrix and then use modular arithmetic (again with 
p = 101) to compute its rank; the result is 84, so the dimension of the nullspace is 
15316. To decide whether any nullspace vector represents a new identity in degree 9, 
we need to determine the dimension of the subspace of the nullspace which consists 
of all consequences of the ternary recombination identity R = _R(a, b, c, d, e, /, g) 
from Theorem 1101 To obtain consequences in degree 9, we can either 
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(i) replace one of the variables x £ {a, b, c, d, e, /, g} by the triple {x, h, i} where 

h and i are two new variables, or 
(ii) embed the entire identity R in a triple {R, h, i}. 

This gives 8 identities in degree 9 which generate the Sg-module of all consequences 
of the identity R in degree 7: 

R({a, h, i}, b, c, d, e, f, g), R(a, {b, h, i}, c, d, e, /, g), 

R(a, b, {c, h, i}, d, e, f, g), R(a, b, c, {d, h, i}, e, /, g), 

R(a, b, c, d, {e, h, i}, f, g), R(a, b, c, d, e, {/, h, i},g), 

R{a, b, c, d, e, f, {g, h, i}), {R(a, b, c, d, e, f, g), h, i}. 

We need to "straighten" the terms of these identities: that is, we use the complete 
symmetry identity to replace each monomial by an equivalent monomial in one of 
the four association types listed above. 

We now apply all 9! permutations of a, b, c, d, e, /, g, h, i to each of the 
8 consequences of R in degree 9, store these permuted identities in the rows of a 
matrix with 15400 columns, and compute the rank of the matrix. For £ = 1, 2, . . . , 8 
the dimension of the subspace generated by the first I consequences of R is as 
follows: 

14071, 15036, 15036, 15036, 15316, 15316, 15316, 15316. 

Since 15316 is also the dimension of the nullspace of the expansion matrix, it fol- 
lows that every polynomial identity in degree 9 satisfied by ternary intermolecular 
recombination is a consequence of the ternary recombination identity in degree 7. 
Hence there are no new identities in degree 9. 

4. Conclusion 

The results of Bremner [T|, Sverchkov [5] and the present paper are consistent 
with the following conjectures for the general case of n-ary intermolecular recom- 
bination. 

Conjecture 1. Every polynomial identity satisfied by n-ary intermolecular recom- 
bination in degree 2n—l (that is, in which each term involves two applications of 
the n-ary operation) is a consequence of complete symmetry in degree n. 

Conjecture 2. Every polynomial identity satisfied by n-ary intermolecular recom- 
bination in degree 3n— 2 (that is, in which each term involves three applications of 
the n-ary operation) is a consequence of complete symmetry in degree n together 
with an identity R in degree 3n— 2 in which every coefficient equals ±1 and which 
implies that every monomial in the second association type, 

{ { (Zl, . . . , a n }, { a n +l, • ■ • , 02n }: 0>2n+l, ■ • ■ , fl3n-2 }, 

equals a linear combination of monomials in the first association type, 

{ { { Ol, . . . , a n }, a n+ i, . . . , a2„-i }, 02«, . . . , fl3n-2 }■ 

Conjecture 3. Every polynomial identity satisfied by n-ary intermolecular recom- 
bination, with no restriction on the degree, is a consequence of complete symmetry 
in degree n together with the identity R from Conjecture^ in degree in— 2. 

Resolving these conjectures is the most important open problem in the theory 
of n-ary intermolecular recombination. 
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Remark 12. Shortly before the final version of this paper was sent to the editors, 
I received Sverchkov's preprint |6], which contains a complete proof of all three 
Conjectures for all n, with an explicit identity R as required by Conjecture^ 
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